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PARTICULAR SOLUTIONS OF THE INVERSE PROBLEM EOR ID 
VLASOV-MAXWELL EQUILIBRIA USING HERMITE POLYNOMIALS 

O. ALLANSON, T. NEUKIRCH, S. TROSCHEIT AND F. WILSON 


Abstract. We present the solution to an inverse problem arising in the context of finding 
a distribution function for a specific collisionless plasma equilibrium. The inverse prob¬ 
lem involves the solution of two integral equations, each having the form of a Weierstrass 
transform. We prove that inverting the Weierstrass transform using Hermite polynomi¬ 
als leads to convergent infinite series. We also comment on the non-negativity of the 
distribution function, with more detail on this in Allanson et al, Journal of Plasma 
Physics, vol. 82 (03), 2016. Whilst applied to a specific magnetic field, the inversion 
techniques used in this paper (as well as the derived convergence criteria and discussion 
of non-negativity) are of a general nature, and are applicable to other smooth pressure 
functions. 


1. Introduction 


Plasma equilibria are of great importance for laboratory, space and astrophysical appli¬ 
cations, see [221131] for example. In this paper we present results arising from an inverse 
problem in the context of collisionless and non-relativistic plasma equilibria, namely that 
of calculating the 1-particle distribution function, given a prescribed magnetic field. We 
consider a system depending only on one Cartesian spatial coordinate {x, y, z), namely the 
z-coordinate. In line with previous work [8l ll8l[26] . we will make several simplifying as¬ 
sumptions. Firstly, we assume a two species plasma of ions and electrons with equal and 
opposite charge of magnitude e, such that the charge density vanishes. Then the electric 
held E vanishes, and Poisson’s equation is automatically satished. A second assumption 
we make is that the magnetic held B only has two non-vanishing components, namely that 
B( 2 :) = By{z),0). For an equilibrium, one then has to solve the coupled Vlasov- 

Maxwell equations 


( 1 ) 

( 2 ) 

( 3 ) 


d? At 


dz^ 

dz^ 


+ —(vx B) -Vv fs = 0, 
oz rris 

= Pojx = j Vx fs Av, 


= Pojy = Po^qs / Vyfsd^V, 


where fs{z,Vx,Vy,Vz) (> 0) is the 1-particle distribution function (DF) for particle species 
s with charge Qs and mass rus, dehned over a AD phase-space of z and v. We have also 
introduced the vector potential A(z) = {A^ (z), Ay (z) , 0) with B = Vx A, thus automatically 
satisfying V • B = 0, and the electric current density j = (jx, Jy,0) = V x B/^o, with its 
dehnition in terms of velocity moments of the distribution function fs given by the x and y 
components of Ampere’s law, and o, respectively. 

The Vlasov equation, o, is a hrst order partial differential equation and the character¬ 
istics are the equations of motion of a charged particle in the magnetic held B. Due to the 
symmetries of the problem, one can immediately hnd three constants of motion: namely 
the Hamiltonian (energy) of a particle of species s, Hg = TOsU^/ 2, and the two canonical 
momenta, Pxs = rrisVx + qsAx and Pys = nisVy + qsAy. Any non-negative and differentiable 
function fs{Hs,Pxs,Pys) with the property that its velocity space moments of all order exist, 
is then an acceptable solution of © , see m- 
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The ‘forward probleni’ of finding Vlasov-Maxwell equilibria of the type just described is 
defined by choosing a distribution function fs[Hs,Pxs,Pys) for each species s and using it 
to determine the current density components jx{Ax,Ay) and jy{Ax, Ay). The final step in 
completing the solution is solving the two coupled ordinary differential equations ([2]) and 

If one would like to determine distribution functions for a given magnetic field B(z), one 
has to solve the ‘inverse problem’ of collisionless plasma equilibrium theory. Channell [5] 
has shown that it is possible to formulate this inverse problem as an integral equation using 
one component of the pressure tensor, which in the case of our choice of coordinates is the 
P 22 -component, defined by 

Pzzi^Ax^ Ay) = ^ ^ 171 s I ^z fsi^s: Pxs: Pys) ^ • 


Pzz is formally found by taking the moment of /s by — Uzs)'^, for Uzs the fluid 

velocity in the z direction. However, as in [8] we will assume that the form of fs as a function 
of its three variables is the same for both particle species and that in particular 


fs{Hs,Pxs,Pys) 


{s/^Vth.s)^ 


exp{-f3sHs) gs{p XS 5 Pys)- 


This is an even function of Vz and hence Uzs = 0, since 


1 

ns{r,t) 


V fs d^v, 


for Ug the number density of species s. The parameter j3s = (.’’^sVth s)~^ normalises the 
Hamiltonian, with Vth,s the thermal velocity of species s. uqs is a species dependent con¬ 
stant and ps is an unknown function that has to be determined by solving the inverse 
problem. This is the main calculation of this paper, ps will differ between species only 
due to its dependence on species dependent parameters such as ms, Vth,s or ps, but not via 
its dependence upon the two canonical momenta Pxs and Pys- Assuming that the function 
PzziAx, Ay) is known, the inverse problem is defined by the integral equation 


Pzz{Ax, Ay) — 


fde + Pi 


^Os 


( 4 ) 


PePi ‘^'n-ms'^vth.s'^ 

/ OO POO 

-OO j —OO 


as explained in [5] , with examples of solutions in [21IH1IIZ1I111I3S] • Charge neutrality requires 
the RHS of (01) to be independent of species, see [11S1II71I35] for examples and discussion. 

In astrophysics, there is a particular interest in force-free magnetic fields, see for example 
[2^ . These are fields such that the electric current density j is everywhere parallel to B. 
This can be described mathematically by j = a(r)B, causing the Lorentz force to vanish, 
j X B = 0. The force-free field is said to be nonlinear if a varies with position, linear 
otherwise and potential if a = 0. This paper is concerned with a particular nonlinear 
force-free configuration; the ‘Force-Free Harris Sheet,’ hereafter referred to as FFHS, see 
[T7 l [T ^ [26ll^[36]. 

B = Ho (tanh(z/L), sech(z/L), 0). 

The FFHS describes a current sheet - a structure of particular interest in reconnection 
studies [55] - with L the ‘width’ of the current sheet and Bq the constant magnitude of the 
magnetic field. A force-free equilibrium in Cartesian geometry implies that Pzz must have 
zero spatial gradient (force balance) [5S]. It is shown in m that the Pzz must also satisfy 
Ampere’s law in the form 

^ dPzz 

dz"^ dA 

As shown in ng, Ampere’s law only allows a force-free solution if the function Pzz{Ax, Ay) 
is constant along this solution, i.e. if the solution is also a isocontour of Pzz over the Ax- 
Aj,-plane. In the nonlinear force-free cases found by m and [5] the function Pzz{Ax, Ay) 
had the form Pzz = Pi{Ax) -I- P 2 {Ay) (called the ‘summative form’ from now on). 
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However, as shown in [18] and discussed in HZ], Ampere’s law admits an infinite number of 
pressure functions for the same force-free equilibrium. Once a Pzz with the correct properties 
has been found one can define another Pzz giving rise to the same current density (for the 
force-free solution) by using the nonlinear transformation 

(5) P ..=^\ Pff )-^ HP ..). 


Here any differentiable, non-constant function ip can be used, such that the RHS is positive, 
with Pf f the pressure, Pzz , evaluated at the force-free vector potential Ay/. 

Obviously, even if the integral equation (|3|) can be solved for the original function Pzz{Ax, Ay) 
it is by no means clear that this is possible for the transformed function Pzz- Usually one 
would expect that solving (|3|) for gs is much more difficult after the transformation to Pzz ■ 
Hence, it is so far unclear whether the transformation property is just a mathematical cu¬ 
riosity or of any practical use. 

In this paper we show for the first time that the transformation property of Pzz can 
actually be useful. As a particularly interesting choice of transformation, this paper deals 
with exponential transformation 


( 6 ) 


Ip (Pzz) = exp 



Pff) 


for Pq a positive constant. It is readily seen that Pzz\Aff = Po and so free manipulation of 
the constant pressure is possible. This is of particular interest because it potentially allows 
us to freely choose the value of the ratio between the pressure and the magnetic pressure 
i?^/(2yto), the so-called plasma beta. This is not possible in the cases in which Pzz is of 
summative form, which inevitably have a lower bound of one for the plasma beta. Note that 
if Pzz is of a ‘summative form’ prior to transformation, then it will be of a ‘multiplicative 
form’ afterwards, since 


exp(Pi(A„;) -H P2{Ay)) = Pl{AjP)P2{Ay). 


The main aim of this paper is to show that there is a mathematically meaningful method of 
solving Q for the exponentiated form of the pressure function used in |17j for the force-free 
Harris sheet. 

This paper is organised as follows. In Section 2 we find that by transforming the summa¬ 
tive pressure function of m according to and dH), we are able to calculate a new DF, 
i.e. we solve the inverse problem defined in (|4|) - with exponentiated functions on the LHS - 
for the unknown function gs ■ We use infinite series of Hermite polynomials in the canonical 
momenta to do this, and prove that the series are convergent and bounded in Section 3. The 
necessary matching between the microscopic parameters of the distribution function and the 
macroscopic parameters of the force-free equilibrium is given in Appendix A. In Appendix 
B we comment on the non-negativity of the distribution function. 


2. Formal solution to the inverse problem 


The ‘summative’ pressure used in m for a FFHS equilibrium is of the form 


(7) 


Pzz{Ax, Ay) 



1 

-cos 


2Ax 


-I- exp 


2Ay 

BoL 


+ Pb. 


Pb > i?g/(4/j,o) is a constant that ensures positivity of Pzz and /ip is the magnetic perme¬ 
ability in vacuo. This is the function that we exponentiate according to (O and (|n|). Clearly 
the resultant function will be ‘multiplicative.’ We choose a pressure function and gs function 
of the form 


Pzz = nos exp - 


gs = exp 


1 /3e + /3i 


WpiJ 


Pl{Ax)P2iAy), 


2Ppi 


9is{p XS )g2s{Pys), 
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such that ppi = 2/4 o-Po/Sq is the plasma beta, a dimensionless constant, and not to be 
confused with the thermodynamic beta /3s. Separation of variables in (|3]) leads to 


(8) 

A (A,) = 

1 

POO 

> —00 

V^rnsVth,s . 

(9) 

A (A,) = 

1 

pOO 

/ —00 

V^rnsVth,s - 


The pressure tensor components are now written as IZ^ integral transforms of the unknown 
parts of the distribution function. The Weierstrass transform, $(a:) of (j){y), is defined by 

(10) $(a;) = yV[4>\{x) = —= / e~*'^~y^''/'^(/){y)dy, 

VdTT J_oo 

see [5]. This is also known as the Gauss transform, Gauss-Weiertrass transform and the 
Hille transform, m- Intimately related to the heat/diffusion equation, $(a;) represents 
the temperature profile of an infinite rod one second after it was (/'(x), [34], implying that 
the Weierstrass transform of a positive function is itself a positive function. Pi and P 2 
are expressed as Weierstrass transforms of gis and ( 72 s in ([S]) and ®, respectively, give 

r.2 

or take some constant factors. Formally, the operator for the inverse transform is e~^ 
with D the differential operator and the exponential suitably interpreted. See [Tl]|35] for 
two different interpretations of this operator. A third and perhaps more computationally 
‘practical’ method employs Hermite polynomials, Hn { x ), see |S], defined here by a Rodrigues 
formula; 

Hnix) = —e- . 

The Weierstrass transform of the Hermite polynomial 7Z„(y/2) is a;". Hence if one knows 
the coefficients of the Maclaurin expansion of $(a;); 

00 

j=0 

then the Weierstrass transform can immediately be inverted to obtain the formal expansion 

00 

(11) (t>iy) 

j=0 

For this method to be useful in our problem, the pressure function defined by ([S]) and 
must have a convergent Maclaurin expansion, with its coefficients allowing the Hermite 
series to converge. For the use of Hermite polynomials in Vlasov equilibrium studies see 
|7l|8j[T4l[T5l|T9l[33] . The question of convergence of a formal solution made up of an infinite 
series of Hermite polynomials was raised in [TO] , but not answered. In [53] , the DF is assumed 
to be a Maxwellian, multiplied by a sum of Hermite polynomials in velocity and some 
arbitrary dependence on the vector potential, to be determined later on. Our approach sheds 
some light on the relationships they derived, whilst beginning from the constants of motion 
approach. Finally, we note that one of the existing nonlinear force-free Vlasov-Maxwell 
equilibria known, is based on an eigenfunction of the Weierstrass transform, m- 
We shall follow a method of the same ilk as Ghannell, that is to expand the ps function 
in Hermite polynomials, with as yet undetermined coefficients. This method is chosen as we 
believe that it provides the neatest route to deriving ‘consistency relations’. These consis¬ 
tency relations - derived in Appendix A - are necessary to link the macroscopic parameters 
that describe properties of the current sheet {Bq, Pq, L) to the microscopic ones that relate 
to the individual particles {rus, Vth,s, Ps, i^os)- 

2.1. The LHS of ([4]): Maclaurin expansion of the transformed pressure There 
is a result from combinatorics due to Eric Temple Bell, that allows one to to extract the 
coefficients of a power series, f{x), that is itself the exponential of a known power series. 
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h{x), see Bell’s original paper [4]. If f{x) and h{x) are defined 
(12) f{x) = e- 


Hx) = —^Cr, 

ml 

m—l 


then we can use ‘Complete Bell polynomials’, also known as ‘Exponential Bell Polynomials’ 
and hereafter referred to as CBP’s, to write f{x) as 


(13) 


1 


f{^) = E ^^-(Cl,C2,...,Cm)x'". 

m \ 


m—0 


^m(Ci 5 C 25 ■•■Cm) is the m^^ CBP. Instructive references on CBP’s can be found in | 

For this paper, the Maclaurin coefficients for the exponential and cosine functions of equation 
0 are used as the arguments of the CBPs. These CBPs are used to form the Maclaurin 
coefficients of Pi and P 2 as in equation (IT51) . Let us now calculate the Maclaurin expansion 
of Pi using CBPs. First expand the cosine: 


Pi = exp 


1 


2Ppi 


cos 


( 2 A,) 


= exp 


1 




exp 


^ ml 


such that a 2 i-i = 0 and 0 : 2 ; = (—1)^2^* ^/Ppi for / G N. is normalised by BqL to give 
Ax- Now that we have an expression of the form of equation m, we can use the CBPs to 
express Pi as a power series in A^: 


Pi = 

such that 

(14) a 2 m = exp 


/ 1 \ 1 . . 

(^) ^ 7^-^Um(ai,a2,-,a2m)i^'" = E 

' m=0 m=0 

^_j^^m22m — 1 


2ppiJ (2m)! 


Ppr'^’Ppi 


There is a standard result, [T^, of CBPs that Fm(aCi, a^C 2 ) a™C™) = a’”Lm(Ci) Cm)- 

Hence we can simplify equation (1141) to give 


(15) 


a 2 m = exp 


1 


(- 1)™2 


mo2m 


-r,, 


0,:rP,0,...,0 ^ 


2ppi) (2m)! V”’2/3pz''’2/3pz 

The odd terms are all zero and hence not included. This makes sense since exp(cos(a;)) is 
an even function. 

Next, we find the Maclaurin expansion of P 2 {Ay) using the same method as above. 


P 2 = exp 


Ppl 


exp 


1 

Ppi 


e-p(ES^; 


\n—l 


such that jn = 2"//3pi. Ay is normalised by BqL to give Ay. Making use of equation (fUl) . 
we can now express P 2 as a power series in Ay\ 


(16) 

such that 


P 2 = exp 


.j \ 00 .j 00 


n—0 


n=0 


bn = exp 


1 \ 1 


Yn 


2" "N 


- . 


^PpiJ n\ ” \l3pi ’ ■■■’ PpiJ 
Noticing a geometric factor, we can simplify the CBPs as before to give 


(17) 


bn = exp 


1 




1 


1 


^PpiJ nl ” \l3pi ’ ■■■’ PpiJ ' 

To conclude, the new pressure function P^z - calculated by exponentiating Pzz as in 0, 
according to ® and ([6]) - has the Maclaurin expansion 

-1 


(18) 


Pzz = nos exp 


2f3piJ Mi 




m—0 


n=0 
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with a 2 m and bn defined in (HID and Cll) respectively. 

2.2. The RHS of Q: Hermite expansion of gg. We make the ansatz that gg is de¬ 
scribed by the following expansions, in line with the preceding discussion on inversion of the 
Weierstrass transform; 



for currently unknown species dependent coefficients C 2 m,s and Dns- We cannot simply 
‘read-off’ the coefficients of expansion as in m, since our integral equations are not quite 
in the ‘perfect form’ of (nni). Upon computing the integrals of equations (|8]) and ([9]) with 
the above expansions for gg, we have 



and hence 



This result appears species-dependent. However, in accordance with [S1III1I2S] we have to 
fix the RHS to be both species-independent, and to match with the pressure function that 
maintains equilibrium with the FFHS, i.e. dT51) . The results of this calculation are in the 
Appendix. 

We are now able to formally write the DF: 



(y/^Vth.s) 



(19) 


with C 2 m,s and Dns given in the Appendix. 


3. Convergence and boundedness of the distribution 


Here we prove that the Hermite series in m converge. A sufficient, general convergence 
criterion on the coefficients of expansion for is also established, that could be applicable 
for other series/equilibria. 

Theorem 1. The infinite sums 



converge V Pxs, Pys S K. and nig, Vth,s S with C 2 m,s and Dns defined in the Appendix. 

We shall first prove the convergence of g 2 s using a combination of the comparison and 
the ratio test, and then prove convergence of gig by comparison with g 2 s'. 
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Proof. Rather than using the CBP formulation of the Maclaurin expansion of the pressure, 
we shall calculate the expansion ‘explicitly’. By twice using the Maclaurin expansion of the 
exponential, we see that 


exp 


1 

Ppi 


.,2A, 


s.t. K = 


n—0 


Using a result from the Appendix we have 

OO 

^2s(Pys) — ^ ^ DnsHn 


Pys 


n—0 


( 20 ) 


= 


n—0 


E 


V^lTlsVth.s 

n / 

H, 


Pys 


\/2i 


'^s'^th.s 


An upper bound on Hermite polynomials is provided by the identity (see e.g. [50] ! 
(21) \Hjix)\ <k^/j'.2^/^exp{x^/2) s.t. /c = 1.086435. 

This upper bound implies 


Pys 


y/2', 




< kbnS'fv n\ exp 


( pIs 


Working on the level of the series composed of upper bounds, the ratio test requires 

bn+l 


( 22 ) 


lim 

n—¥oo 


\/n+l = lim \/n + 1 < l/5s 

n^oo Ofi 


for a given Ss S K+. Note that at this stage, we are working with a general bn, and hence 
equation (E21) is a sufficient criterion for the convergence of a Hermite polynomial series of 
the type in equation (l20ll . corresponding to a pressure expansion of the type of equation 
(US- The next step is to analyse the asymptotic behaviour of Explicit expansion 

of the bn coefficients gives 


Wra+l 


/bn = 


OO 


2Jjr. 




/ 2 


2(^n 


n + 1 


0!/3pi 


+ 


2"^2 


, OO 

E 

3 

2^3 




1!E 


2!E 


2‘on 


V l!/3pi 

(1 

Ppi 


+ 


2^2 


Son 


2'‘*3 


21/32 


pi 
2r\n 


31/3.3 


n + 1 


^ + 2 


2^2 


2!/3^z 


+ 3 


pi 


2on 


2^2 


V ll/3pz 21/3.2 


pi 


233” 


The kth ‘partial sum’ of this fraction has the form 

Pi + 2p2 + 3p3 4 


Tfc = 


kpk 


Pl+P2+P?. + ■■■ 


\ 

) 

\ 

J 


with Pi X 1/il, where we write g h to mean g/h and h/g are bounded away from 0. Now 
since the denominator of the pi increase super-exponentially (factorially) we have ipi x pi 
and hence 

OO OO 

0 < ipi < OO and 0 < pi < oo, 

i—1 i—1 

giving Tk —>■ Too G R’*’ and, more specifically, Too x 1 in n. Therefore bn+i/bn = Too/in+l) x 
1/n. That is to say bn+i/bn grows asymptotically like 1/n. We then get 

lim -v/n + 1 = lim 1/y/n < 1/Ss V<5sGK^. 

n—^oo Oji n—¥oo 
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Therefore a series with each term at least as large as those of g 2 s{Pys) converges with an 
infinite radius of convergence (Vds) by the ratio test. Hence g 2 s{Pys) also converges V 5s and 
Pxs by the comparison test. We shall now prove convergence of gis, by comparison with g 2 s- 
By explicitly using the Maclaurin expansion of the exponential, and then the power series 
representation for cos” x from m- 


COS^” X = 


COS^" ^ X = 


22r, 


22 


- n—1 


fc =0 


2n- 1 
k 


cos(2(n — k)x) + 

cos((2n — 2k— l)x), 


one can calculate 


( 1 \ oo 

— cos (2As;^ ) ^ E 

' m—0 

The zeroth coefficient is given by qq = exp (l/(2/3p;)), and the rest are 


^2m — 


2(-ir 


■EE 


(j - 


for Jk = {2k + 1, 2fc + 2,...} and m 0. By rearranging the order of summation, a 2 m can 
be written 

L(j-i)/2j 


a-2m — 


■E 


1 


2(-l)^ 

( 2 to )! 


E 


(j - 2fc)^™, 


j — l " - • k—0 

where [xj is the floor function, denoting the greatest integer less than or equal to x. Using 
a result from the Appendix, we have 


9i 


\{,Pxs) — ^ ^ C27n,s^2va 


Pxs 


m—0 


— ^ ^ 0'2m 

m—0 


E 

V2 


y\/2msVth,s 

2m / 

H2m 


Pxs 


V2: 


'msVth^s 


Recognising an upper bound in the expression for a 2 m] 


l(i-i)/2j ... j 

n—0 ^ ^ n=0 


gives 


2(_l)m ^ 2J + lj2™ 
a2m ^ 77. 7i~ / . 


= 2 


j2m 


(2m)! ^^j\V{2f3piy (2m)! ^ j!(2/3p0 


■E 


< 


2 j2m 


■E 


(2m)!2^j!(2/3p0^’ 

■ 2^ -^ ^2m 


(2m)! ^ jiy 


Hence we now have an upper bound on a 2 m for m ^ 0 and we know that a 2 m+i = 0, and 
so is bounded above by b 2 m+i- Note also that uq < &o- Hence, each term in our series for 
gisiPxs) is bounded above by a series known to converge Vds according to 

(-^) ^(:3) 

So by the comparison test, we can now say that gis (Pxs) is a convergent series. This 
completes the proof of Theorem [T] □ 
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Note that it is not sufficient for the distribution to be merely convergent. It must also be 
bounded by a finite constant over all momentum space. Our DF is bounded in momentum 
space by 



where Sis and S 2 s are finite constants, since gis and g 2 s are convergent. The additional 
exponential factors come from the upper bounds on Hermite polynomials used above m- 
This clearly goes to zero for sufficiently large \pxs\, \Pvs\ and is without singularity. We 
conclude that the distribution is bounded/normalisable. 


4. Summary 


Starting from a one-dimensional, nonlinear force-free current sheet equilibrium, we have 
found a new solution to the inverse problem of calculating non-relativistic, neutral colli¬ 
sionless equilibrium distribution functions for both ions and electrons. This was achieved 
by transforming a known pressure function by exponentiation, exploiting the fact that an 
equilibrium pressure function is non-unique. It was noted that by using ‘Channell’s method’, 
the resulting integral equations expressed the pressure as Weierstrass transforms of the un¬ 
known parts of the distribution. These were inverted by using a known technique, namely 
a matched expansion in Hermite polynomials. Complete Bell polynomials were also used 
to determine the form of the pressure tensor. However, this solution is purely formal, until 
convergence, boundedness and non-negativity are proven. In Section 3, convergence of gs 
was proven by using the ratio and comparison tests, and boundedness was demonstrated. 
For a brief discussion of the non-negativity of the DF, see Appendix B, and for the full 
details see [3]- 

A significant technical problem has also been overcome, since previous force-free solutions 
[1II71I31I3S] could not describe low-beta configurations, which are prevalent in space and 
astrophysical plasmas. Whereas the distribution in this paper can describe such plasmas. 
This could be significant for modelling processes such as reconnection/tearing mode in low- 
beta environments like the Solar Corona or the Scrape-Off Layer in Tokamaks, see [sna 
respectively. Further work in stability analysis and general properties of the DF is a priority. 

The spirit of the methods used in this paper are rather general, and could be applied 
to other magnetic fields with smooth pressure functions P{Ax, Ay), force-free or otherwise: 
namely the ‘matching’ of Maclaurin and Hermite expansion coefficients; the convergence 
criteria and the positivity argument. It seems a legitimate project to formally generalise the 
procedures used in this paper, perhaps extending to include a non-zero electric potential. 
This could result in a reliable analytic algorithm for calculating valid solutions to the inverse 
problem for collisionless plasma equilibria, for arbitrary ID, 2-component magnetic fields. 
A similarity between solutions of the inverse problem discussed herein, and Green’s function 
solutions of the Diffusion Equation has also been established. 
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Appendix A. Consistency of the micro and macroscopic descriptions 

The result of the Maclaurin expansion is to be able to write the pressure tensor in terms 
of macroscopic parameters as 


m=0 n—0 

To ensure neutrality ni{Aj:,Ay) = UeiA^, Ay), we require the RHS of equation (|4]) to be 
identical for ions and electrons, see min]. We can let uqs be independent of species, i.e. 
noi = noe = no- To fix the other parameters of the DF, we insist that the Maclaurin 
expansion of the pressure tensor above matches up with those found in Section [2.21 
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and Pq = is the species-dependent magnetisation parameter, [13], also used in 

gyrokinetic theory as the fundamental ordering parameter, |20| 
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It is the ratio of the Larmor radius, p, = to the characteristic length scale of the 

system, L. Q, = gsBo/m, is the gyrofrequency of particle species s and T, is the temper¬ 
ature, ksT, = rUsV^f, When S, << 1 then particle species s is highly magnetised and a 
guiding centre approximation will be applicable for that species, see m- The conditions 
derived in this Appendix are critical for ‘fixing’ the DF, and making a link between the 
macroscopic description of the current sheet, with the microscopic one of particles. 

Appendix B. Non-negativity of the Distribution 

Since a DF represents a probability density, then it clearly must be non-negative over all 
phase space. The non-negativity of DFs that are partially represented by sums over Hermite 
Polynomials has been raised by previous authors, see for example wm- We refer the reader 
to a recently published article, [3], for further discussion of non-negativity of g, functions, 
represented by (possibly) infinite sums over Hermite Polynomials. To paraphrase the results 
therein, it is claimed that - provided g, is differentiable and convergent - the g, function 
will be positive for a sufficiently magnetised plasma, i.e. for 0 < dg < Sc < cxj, for some 
critical value of the magnetisation parameter, Sc- 
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